Ballistic- Transport-Induced Negative Differential Thermal Resistance 
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Abstract 

Using nonequilibrium molecular dynamics simulations, we study the temperature dependence of the negative differential 
thermal resistance that appears in two-segment Frenkel-Kontorova lattices. We apply the theoretical method based on Landauer 
equation to obtain the relationship between the heat current and the temperature, which states a fundamental interpret about 
the underlying physical mechanism of the negative differential thermal resistance. The temperature profiles and transport 
coefficients are demonstrated to explain the crossover from diffusive to ballistic transport. The finite size effect is also discussed. 
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Negative differential thermal resistance (NDTR) is a 
property of certain materials connected to two heat baths 
in which, over certain temperature ranges, heat current 
is a decreasing function of the temperature difference be- 
tween the two heat baths, i.e., dJ/dT < 0, here J and T 
are the heat current and the temperature, respectively. 
Li et al had investigated different kinds of systems to re- 
veal that negative differential thermal resistance is not 
a happenstance and indeed take place in low tempera- 
ture nonlinear lattices for a certain range of parameters 
H il il- It is naturally expected that NDTR ef- 
fects may lead to an impressive technological innovation 
and even the appearance of the thermal transistor, which 
maybe change our world thoroughly in future, just like 
electronic transistor and other relevant devices have done 
in the past half-century 1] [5] [6]. 

Although thermal transistor is confirmed by nonequi- 
librium molecular dynamics simulations, the practicable 
counterpart using fabricated materials is not produced 
up to now [ij. As the main physical mechanism of ther- 
mal transistor, negative differential thermal resistance 
(NDTR) effect as well as asymmetric heat conductance 
still need to be understood deeply The phenomenon 
of NDTR is firstly understood from the mismatch be- 
tween the phonon bands of the two interface particles 
In Ref. 0] and [oj, the authors reported that the 
phonon bands are independence of the system size but 
the NDTR disappears for the large system. Segal has 
suggested that NDTR shows up when the molecular is 
strongly coupled to the thermal baths in a asymmetric 
system However, the NDTR in the absence of asym- 
metry is still not clear Therefore, the mismatch of 
the phonon bands and asymmetry are not a real physi- 
cal mechanism of the NDTR. To reveal the real physical 
mechanism of the NDTR is still one of new and challeng- 



ing problems about thermal transport. 

Macroscopic Fourier's law, J = — trVT, that connects 
heat current with thermal transport on the microscopic 
scale, is an empirical law based on observation, where VT 
is the temperature gradient and a is the thermal conduc- 
tivity. It states that the heat current through a material 
is proportional to the negative temperature gradient. On 
the other hand, according to microscopic thermodynam- 
ics, heat current at low temperature is proportional to the 
number of phonons [11]. A phonon is a quantized mode 
of vibration occurring in a rig id crystal lattice, such as 



the atomic lattice of a solid . The Bose-Einstein prob- 



ability distribution for phonons, based on statistical me- 
chanics concepts for thermal equilibrium, determines the 
number of phonons. Thus, the heat current is also deter- 
mined by the Bose-Einstein probability distribution. In 
this paper, applying an analytical method from Landauer 
equation as well as nonequilibrium molecular dynamics 
simulations, we reveal the dependence of the NDTR on 
the temperature and the crossover of two thermal trans- 
port processes: the ballistic transport and the diffusive 
transport We will also investigate the system size 

dependence of the NDTR. The temperature profile and 
the transport coefficient will be calculated to characterize 
the ballistic transport and the diffusive transport. 

The nonlinear lattices that we use in this letter consist 
of two segments, left segment (L) and right segment (R). 
Each segment is a Frenkel-Kontorova (F-K) lattice. Seg- 
ment L and R are coupled via a spring of constant ifmt ■ 
The total Hamiltonian of the model is 



H — Hl + Hb + Hi, 



(1) 
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and the Hamiltonian of each segment can be written as 
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with qM,i and pM,i denote the displacement from equi- 
hbrium position and the conjugate momentum of the 
particle in segment M, where M stands for L or 
R. The parameters K and V are the harmonic spring 
constant and the strength of the external potential of 
the FK lattice, respectively. We couple the last parti- 
cle of segment L and R via a harmonic spring. Thus, 
Hint = {qL,N - QR.Nf . We set m a = 1, Kl = 
1.0, Kr = 0.2, Vl = 5.0, Vr = 1.0, K,„t = 0.05. 

In our simulations we use fixed boundary condition and 
the chain is connected to two heat baths at temperature 
Tl and Tr. We use the Nose-Hoover heat baths and in- 
tegrate the equations of motion by using the fourth-order 
Runge-Kutta algorithm [3] The local temper- 

ature is defined as Ti = (pf)- The local heat flux is de- 
fined as ji = KM{pi{qi — Qi-i)), and the total heat flux is 
J = Nj. The simulations are performed long enough to 
allow the system to reach a steady state in which the lo- 
cal heat flux is constant along the chain. For the sake of 
comparison, we define a heat current ratio, Jr = J/ Jmax, 
in which Jmax is the maximum heat current under a fixed 
temperature Tr of the right heat bath. The transport co- 
efficient is an important quantity for characteri zing the 
transport mode of a thermal transport process [l7| [isj 
[l^. The thermal conductance evaluated as cr = Nj/ AT 
represents an effective transport coefficient that includes 
both boundary and bulk resistances [20t . 

Figure 1 displays a typical negative differential thermal 
resistance [l[ . When the temperature difference increases 
(i.e., Tl {< Tr) decreases), the heat current increases 
firstly and then decreases. The former is positive differ- 
ential thermal resistance (PDTR) and the latter negative 
differential thermal resistance. It would be much inter- 
esting to show the thermal conductance dependent be- 
havior of NDTR, which is also presented in Fig.l. It is 
clearly that the NDTR corresponds to the ballistic regime 
and PDTR the diffusive regime. 

The phenomenon of NDTR can be understood from 
the theoretical approach. We use Landauer equation 
to calculate the heat flux. For example, in the ballis- 
tic regime the phonon heat flux can be also calculated 
through a Landauer type expression [2l|: 



J= dLUUJx{uj)[nR{uj) ~ TiLiuj)], 



(3) 



where x(a;) is the temperature independent transmission 
coefficient for phonons of frequency lu. Here UMii-^) — 
(^e^M^ _ i)-i^ fj^ = l/ksTM, (M = L,i?) is the Bose- 
Einstein distribution characterizing the heat baths. 

As reported in Ref. ^ , under the assumption of weak 
system-bath interactions and when going into the Marko- 



FIG. 1: The temperature dependence of heat current ratio 
(blank squares) and the corresponding temperature depen- 
dence of the thermal conductance (solid circles). The sys- 
tem parameters are Nl = Nr ~ 50, Vl = 5.0, Vr — 1.0, 
Kl = 1.0, Kr = 0.2, K,„t = 0.05, and Tr = 0.21. 



vian limit, the probabilities P„ to occupy the n state of 
the phonon are found to satisfy the master equation 

Pn = {n+l)kdPn+i+nkuPn-i-[nkd + {n+l)ku]Pn, (4) 

here the occupations are normalized ^ P„ — 1, and k^ 
and ku are the vibrational relaxation and excitation rates, 
respectively. 

When going into weak system-bath interaction, kd and 
ku can satisfy: 

kd--kL + kR, fc„ = A:Le-^^"° + kRe-"'^'^" , (5) 

with fcjvf = TM{^^)[^ + nM{uJo)], where nM{i^) = {ef"-''^ — 

here, m, luq, a^are molecular oscillator mass, frequency 
and coupling between the system and the heat baths, 
respectively. 

The thermal properties of our model are obtained from 
the stationary state solution of Eq. (4) . The thermal flux 
is given by 

J = LUoY.nikLPn-kLPn-ie-^^'^°). (6) 

here positive sign denotes current flowing from right to 
left. In this equation the first term indicates the thermal 
flux flowing from the L chains into the L heat bath. The 
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second term gives the oppositely flowing flux from L heat 
bath to the chain. The thermal flux could be equivalently 
calculated at the R chain. 

So for the asymmetry system, which is similar to the 
case of a highly anharmonic molecule coupled-possible 
asymmetrically-but linearly, to two heat baths of differ- 
ent temperatures. Here, we simulate strong anharmonic- 
ity by modeling the anharmonic two-segment chains by a 
two levels system that notes a highly anharmonic vibra- 
tional mode [ij. The Hamiltonian for this model is the 
same as presented in Eq.(l), except that we take ti = 0, 1 
only. Following Eq.(3)- (6) and going into the classical 
limit, the heat current reduces into the simple form 



J 



in which 
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(7) 

+ (8) 
(9) 



We can now clearly seek the main factors which in- 
fluence the heat current. The thermal flux is given 
by multiplying four terms: (1) A symmetric prefactor 
^l^r/O^l + Tfl)- (2) The characteristic frequency and 
effective temperature loq/Ts. (3) The temperature dif- 
ference (Tr, — Tl)- (4) The molecular occupation fac- 
tor [exp(a;o/Ts) + 1]^^ 3- This expression denotes that 
the heat current is mainly a competitive effect between 
the temperature difference and the molecular occupation 
factor. In the case of low temperature (i.e., in the bal- 
listic regime), the molecular occupation factor changes 
obviously with Tg by going into an exponential function, 
while in the case of high temperature (Ts>>a;o, indicates 
the diffusive regime) the molecular occupation factor is 
a constant 1/2. 

Figure 2 shows the simulated results of the tempera- 
ture dependence of NDTR. When Tr is small, there exists 
NDTR. As Tr increases, however, the NDTR disappears. 
As also shown in Fig. 3, a theoretical estimation Eq.(8) 
based on Landauer equation conflrms this temperature 
dependence of NDTR simultaneously. 

Why does the NDTR appear at low temperature and 
disappear at high temperature? This phenomenon can 
be also understood in detail from the crossover from the 
diffusive to ballistic transport. In the case of low tem- 
perature Tfl, as shown in Fig. 4a, when changes from 
Tr to zero, the temperature profiles of the lattices illus- 
trate a zero temperature gradient. This diffusive-ballistic 
transition with temperature difference is also shown in 
Fig. 4b. When the temperature difference is small, the 
thermal conductance does not change with T^. How- 
ever, when the temperature difference is large, the ther- 
mal conductance decreases linearly with Tl decreasing. 



FIG. 2: Heat current ratio as a function of the temperature 
Tl for rR=0.11, 0.31, 0.51, and 1.01. The system parameters 
are Vl = 5.0, Vr = 1.0, Kl = 1.0, Kr = 0.2, Kint = 0.05, 
and Nl=Nr^^O. 



FIG. 3: The temperature dependence of heat current. Results 
are obtained from Eq.(8) with Tfl=200K (fuU), lOOK (dash), 
50K (dot). The remain parameters are ojo = 150meV, Tl = 
l.OmeV, r_R = lAmeV. 
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FIG. 4: (a) Temperature profiles of the lattices for differ- 
ent heat baths for rL=0.02. (b) Thermal conductance as a 
function of the temperature Tl /Tr for different temperatures 
Tr^O.II, 0.21, 0.41, 0.81 and 1.01 (from bottom to top). Re- 
main parameters are the same as for Fig. 2. 

It can be interpreted that the difFusive-ballistic transi- 
tion induces the NDTR, which takes place just when 
the ballistic transport prevails over the diffusive trans- 
port. In the case of high temperature Tr, even when 
Tl is very small, the temperature gradient are non-zero 
and the thermal conductance does not change with T^. 
There exists no diffusive-ballistic transition. Therefore, 
the NDTR does not occur in the case of high temperature 
Tr. 

We would also like to discuss the finite size effect of 
NDTR. As shown in Ref. When the system size 

increases, the phenomenon of NDTR disappears. It is 
interesting to present the system size dependence of ther- 
mal conductance, which is shown in Fig. 5. In the ballis- 
tic regime (i.e., when the system size is much smaller 
than the phonon mean free path), thermal conductance 
increases linearly with the system size. In the diffusive 
regime (i.e., when the system size increases to far larger 
than the phonon mean free path), thermal conductance 
will be independence of the system size. There exists 
a crossover from ballistic to diffusive transport with the 
increase of the system size. When the system goes to 
completely diffusive transport regime, the temperature 
difference does not change the transport mode and then 
the NDTR disappears. 



FIG. 5: Thermal conductance versus lattice length N(— Nl + 
Nr) for Tr = 0.21 and Tl = 0.05 (Triangles), 0.07(Circles), 
and 0.15(Squares). Remain parameters are the same as for 
Fig.2. 

In conclusions, we have found the physical mechanism 
of the negative differential thermal resistance in two- 
segment asymmetric F-K chains through both theoret- 
ical analyses and numerical simulations. From the an- 
alytical relationship between the heat current and the 
temperature, the NDTR effect can be understood from a 
competition between the temperature difference and the 
molecular occupation factor. In the ballistic regime the 
molecular occupation factor mainly depends seriously on 
the temperature. However, in the diffusive regime, the 
molecular occupation factor is a constant. The NDTR 
effect occurs just because the decrease of the molecular 
occupation factor is faster than the increase of the tem- 
perature difference. The main factors that influence the 
transport modes of the phonons are the temperature and 
the system size. We propose that different system struc- 
tures can affect the temperature dependence of NDTR 
through changing the transmission mode of the phonons. 
When there exists a crossover from diffusive to ballis- 
tic transport, the NDTR can be observed, otherwise it 
will not. Due to the rigorous condition of validity: low 
temperature and small system, we still have a long way 
to go before the NDTR effect is produced in fabricated 
materials. 

We would like to thank members of the Centre for Non- 
linear Studies for useful discussions. This work was sup- 
ported in part by grants from the Hong Kong Research 
Grants Council (RGC) and the Hong Kong Baptist Uni- 
versity Faculty Research Grant (FRG). 



4 



[11 

* Electronic address: |wrzhong@hkbu.edu.Iik| j]^2 
[1] L. Wang, and B. W. Li, Phys. Rev. Lett. 99, 177208 

(2007); B. Li, Lei Wang, and Giulio Casati, Appl. Phys. [13 

Lett. 88, 143501 (2006). [14 
[2] W. C. Lo, L. Wang, and B. Li, J. Phys. Soc. JPN 77, 

054402 (2008). [15 
[3] B. Li, L. Wang, and G. Casati, Phys. Rev. Lett. 93, 

184301 (2004). [16 
[4] D. Segal, Phys. Rev. B 73, 205415 (2006). 
[5] J. Bardeen and W. H. Brattain, Phys. Rev. 74, 230 

(1948). [17 
[6] D. Segal, Phys. Rev. E 77, 021103 (2008). [18 
[7] C. W. Chang, D. Okawa, A.Majumdar, and A. Zettl, 

Science 314, 1121 (2006). [I9 
[8] B. Hu, L. Yang, and Y. Zhang, Phys. Rev. Lett. 97, 

124302 (2006). [20 
[9] Z. G. Shao, L. Yang, and B. Hu, Transition From the [21 

Negative Differential Thermal Resistance to the Ordinary 

State, ReportNo. CNS-08-20, 2008. 
[10] B. Q. Ai, Negative differential thermal resistance in sys- 



mmetric F-K single-chains, ReportNo. CNS-10-3, 2008. 
J. A. Reissland, The Physics of Phonons (John Wiley & 
Sons Ltd., London, 1973). 

J. R. Christman, Fundamentals of Solid State Physics 

(John Wiley & Sons Ltd., New York, 1988). 

D. Roy, Phys. Rev. E 77, 062102 (2008). 

S. Nose, J. Chem. Phys. 81, 511 (1984); W. G. Hoover, 

Phys. Rev.A 31, 1695 (1985). 

D. C. Rapaport, The Art of Molecular Dynamics Simu- 
lations, (Cambridge University Press, Cambridge, 2001). 
W. H. Press, S. A. Teukolsky, W. T. Vetterling, and B. 
P. Flannery, Numerical Recipes (Cambridge University 
Press, Cambridge, 1992). 

J. S. Wang, Phys. Rev. Lett. 99, 160601 (2007). 

Y. Xu, J. S. W ang, W. Duan, B. L. Gu, and B. Li, 

larXiv:080 7.3"8T9Vl (2008). 

D. Li, Y. Wu, P. Kim, L. Shi, P. Yang, and A. Majum- 

dara, Appl. Phys. Lett. 83, 2934 (2003). 

S. Lepri, R. Livi, and A. Politi, Phys. Rep. 377, 1 (2003). 

L. G. C. Rego, G. Kirczenow, Phys. Rev. Lett 81, 232 

(1998). 



5 



0.1 



0.01 



1E-3 



1E-4 



Ballistic 
NDTR 



/ 

□ 



dJ/dT<0 / 

□ 



□ 



/ 



1E-5 
0.00 



• •• 



\ 



Diffusive 
PDTR 
dJ/dT>0 



\ 



\ 



0.05 



0.10 



0.15 



0.20 



-4 



2 



-0 



1.0 
0.8 
0.6 
0.4 
0.2 
0.0 



1 

0.1 
0.01 
1E-3 

1E-4 
1E-5 
1E-6 



TR 


0.21 I 

0.81 1 

1.01 j 










(a) 



20 40 60 80 100 



i A i A 4 * * * * * * * ^ 



I r 



/\ 



TR 
—■—0.11 
-•— 0.21 
-A- 0.41 
0.81 
—♦-1.01 



/ 



(b) 



0.0 0.2 0.4 0.6 0.8 1.0 

J jj. 



0.1 











^1 


r* 










/ 






















• / 






■ 






/ ▲ 






























/ 
• 


























/ 




Ballistic 




Diffusive 








(NDTR) 




(PDTR) 

1 



10 100 1000 10000 



N(N+N^ 



